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^ ■ Abstract 

S:^ I We determine the next-to-leading logarithmic (NLL) QCD corrections to the cross section 

^^ [ a{e^e~ — > tiH) for center-of-mass (cm.) energies up to 500 GeV. The dynamics is dominated 

f^ . by nonrelativistic effects, and the summation of terms singular in the relative tt velocity is manda- 

O,' tory to all orders in the strong coupling constant Us using an effective theory. The summations 

(— I I lead to an enhancement of the tree level predictions by about a factor of two and are important 

^ \ for the determination of the top Yukawa coupling. We also study the impact of polarization of 

^ . the electron-positron beams and provide a fast approximation formula for the known 0{as) QCD 

C^ ' fixed-order prediction. 
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I. INTRODUCTION 

The discovery and exploration of the mechanism of mass generation and electroweak 
symmetry breaking is one of the most important tasks of future colhder experiments. Within 
the Standard Model of elementary particle physics (SM) electroweak symmetry breaking 
is realized by the Higgs mechanism which postulates the existence of an electric neutral 
elementary scalar field that interacts with all SM particles carrying nonzero hypercharge and 
weak isospin. Through self-interactions this Higgs field acquires a vacuum expectation value 
V = {V^GfY^^ ^ 246 GeV, Gp being the Fermi constant, which breaks the SU(2)LxU(l)y 
symmetry at high energies down to the electric U(l)em below the symmetry breaking scale 
and leads to nonzero masses of the elementary particles. The Higgs mechanism also predicts 
that the Higgs fields can be produced as a massive Bose particle in collider experiments 
if sufficient energy is provided in the process. The mass of the Higgs boson is expected 
to lie between the current experimental lower limit of 114.4 GeV yj and about 1 TeV. 
Current analyses of electroweak precision observables yield a 95% CL upper indirect bound 
of 186 GeV for the Higgs boson mass |2]. While a Higgs boson with a mass up to 1 TeV can 
be found at the LHC, precise and model-independent measurements of quantum numbers 
and couplings are likely to be restricted to a future e~^e~ Linear Collider J3|, y, |5| such as 
the International Linear Collider (ILC) project. 

The Higgs mechanism predicts that the quark masses mq are related to the quark-Higgs 
Yukawa coupling A^ through the relation mq = \qV. This makes the measurement of the 
Yukawa coupling to the top quark {rrit = 172.5 ± 2.3 GeV J6[) particularly important since 
it is expected to have a high precision. At a future e~^e~ Linear Collider the top Yukawa 
coupling can be measured from the process e'^e~' -^ tiH since the amplitudes describing 
Higgs radiation off the ti pair dominate the cross section. ^ 

For the second phase of the ILC project with cm. energies between 500 GeV and 1 TeV 
and assuming a Higgs mass of around 120 GeV the total cross section a{e^e~ -^ tiH) is at 
the level of 1 — 2 fb and measurements of Xt with experimental errors of around five percent 
are expected [7|, |8[ . The precision motivates the computation of radiative corrections. In the 
approximation that the top quark and the Higgs boson are stable particles^ the tree level 
cross section was determined already some time ago in Refs. 10|. The full set of one- loop 



QCD corrections were obtained in Ref. ll|. Earlier studies using approximations were given 



in Refs. 12, ll3| . One- loop electroweak corrections were studied in Refs. [ij, |l5| and also in 
Ref. ET 



The phase space region where the Higgs energy is close to its upper endpoint, 

EH-E'}, = {s + ml- 4m2)/(2v^) , (1) 



^/s being the center of mass energy, was studied in detail in Ref. |l2|. In the large Higgs 
energy endpoint region the ti pair is forced to become coUinear and to move opposite to the 
Higgs direction in order to maximize the momentum necessary to balance the large Higgs 
momentum, see Fig. QJ Thus the tt invariant mass is close to 2'mt. In this kinematic regime 



^ An indirect measurement through virtual Higgs effects mirfit be also possible at the tt threshold if the 

Higgs mass is close to the present lower experimental limit |3- 
^ For a light Higgs boson this is an excellent approximation. For mn = 115(150) GeV one finds Th = 

0.003(0.017) GeV H 
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FIG. 1: Typical constellation of momenta for the process e^e~ — > tiH in the large Higgs energy 
endpoint region. 



the ti pair is nonrelativistic in its cm. frame and fixed-order QCD perturbation theory in 
powers of Ug leads to singular terms proportional to (ag/v)"' and (ttslnu)"' which have to 
be summed to all order. Here, f = (1 — Amf/Q^Y'^ is the top quark relative velocity 
in the ti cm. frame and Q is the ti invariant mass. In Ref. [l7| these singularities were 
summed at NLL order in a simultaneous expansion in ag and v and also accounting for the 
finite top quark width. The computations were carried out using a nonrelativistic effective 



theory |l9l.l20ll2l| originally developed for the threshold region in the process e^e~ — > ti. Due 
to the large top quark width, Tf ^ 1.5 GeV, the nonrelativistic tt dynamics is protected from 
nonperturbative effects and the summations can be carried out with perturbative methods. 
It was shown in Ref. [TjI that the summation of the singular terms leads to an enhancement of 
the total cross section that needs to be accounted for up to cm. energies of about 700 GeV. 
The impact of the summations increases with the fraction of the phase space where the 
cm. top velocity v is nonrelativistic, i.e. it increases with the Higgs and top quark masses and 
decreases with the cm. energy. A convenient measure for the impact of the nonrelativistic 
summations on the total cross section is the maximal relative velocity of the ti pair which 
is achieved at the low Higgs energy endpoint E^ = ttih, 



For small ti™^^ the summations have a large effect since the available phase space is predom- 
inantly nonrelativistic 



As was already demonstrated in Ref. [l7[, the fixed-order QCD predictions 
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become unreliable for cm. energies up to 500 GeV, which corresponds to the energy available 
during the first phase of the ILC project. For mn = (120, 130, 140) GeV, m^ = 175 GeV, 
and y/s = 500 GeV one has v""^^ = (0.39, 0.32, 0.23) and consequently the entire phase space 
is governed by the nonrelativistic QCD dynamics. The nonrelativistic expansion based on 
the parametric counting «« ~ w -C 1 has to be employed rather then the ag expansion to 
make reliable theoretical predictions for the cross section. Another consequence of small 
^max jg ^\^Q^ ^\^Q cross scctiou for cm. energies up to 500 GeV can be substantially smaller 



than 1 fb due to phase space suppression, which severely restricts statistics. Since the 
singularities proportional to {us/v)^ and Us^^v are large in this case only predictions where 
the nonrelativistic summations are accounted for allow for a realistic assessment of Yukawa 
coupling measurements during the first phase of the ILC project 22, |23| . 



In this work we give a detailed analysis of the total cross section and the Higgs energy 
distribution for the process e+e~ -^ tiH for cm. energies up to 500 GeV accounting for 
QCD effects at NLL order in the nonrelativistic expansion. The approach of Ref. [17| 
developed for descriptions of the large Higgs energy endpoint region is extended to the case 
where the entire phase space is nonrelativistic. We show that our NLL order predictions 
are substantially larger than the known tree level predictions, which have in fact been used 
for experimental simulations studies at 500 GeV in the past 2j|. We also account for 



the possibility of electron-positron beam polarization which can further enhance the cross 
section. Our results significantly affect the prospects for top Yukawa coupling measurements 
during the first phase of the ILC project. 

The content of this paper is organized as follows: In Sec. Owe review the ingredients of 
the factorization formula derived in Ref. jl7[ in the large Higgs energy endpoint region valid 
for large cm. energies. We extend the presentation by also accounting for electron-positron 
beam polarization and by giving a more detailed discussion of the ti final state in the helicity 
basis. In Sec. II I II we discuss the modifications that need to be applied to the factorization 
formula for the case where the full phase space is nonrelativistic. In Sec. II VI we analyze our 
results numerically and Sec. |3 contains the conclusion. 



II. THE LARGE HIGGS ENERGY ENDPOINT REGION 

In the large Higgs energy region Eh ~ E^^ the Higgs energy distribution can be factorized 
into a hard part describing the production of the ti pair and the Higgs boson and in a low- 
energy part describing the nonrelativistic dynamical QCD effects of the ti subsystem. The 
latter are responsible for the singularities proportional to powers of Og/v and a^ In v. The 
factorization formula, valid at NLL order for unpolarized electron-positron beams and top 
quarks, was derived in Ref. [l7[. Accounting for electron-positron beam polarization and 
polarized top quarks the factorization formula for fully polarized electrons and positrons has 
the form 

da .^ .A^ 8iV, [(l+x^-4xi)2-4a;^]i/2 



(;!;'— 2 



H) ) ^3/2 ^2 



i=-l,0,+l 



with 



r/ij _ m'jf _ m^ 



xt = ^, xh = —^, xz = ^. (4) 

s s s 

Here, cq and C(i,i) are the hard singlet and triplet QCD Wilson coefficients which depend on 
the effective theory renormalization parameter u, Fqj_ and F2 .-. _^ are the hard electroweak 
tree-level matching conditions, and C^ is the Green's function of the NLL Schrodinger equa- 
tion of the effective theory for the top quarks. A detailed discussion of these quantities will 
follow shortly. 



The index denotes the hehcity of the electrons, i.e. "— " refers to right-handed positrons 
and left-handed electrons and the index "-I-" refers to left-handed positrons and right-handed 
electrons. Since the electron mass is neglected, the cross section vanishes if both electron 
and positron have the same helicity. For arbitrary polarization P+ of the positrons and P_ 
of the electrons the spectrum reads 
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-(1-P_)(1 + Ph_ 
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(5) 



where the polarization of a beam with Nj^ right-handed particles and N_ left-handed parti- 
cles is defined as 



N+ + N_ 



N,, 



and can take on values between —1 and 1. 

The first two terms is Eq. Q are the hard factors and the third term is the imaginary 
part of the zero-distance Green function of the NLL Schrodinger equation that can be 
derived from the effective theory Lagrangian. The Green function describes the effects of 
the low-energy nonrelativistic dynamics on the ti production rate for the top pair being in 
an S-wave state and does not depend on the polarization of the electron-positron beams. 
It depends on the effective theory renormalization scaling parameter z/ and is proportional 
to the time-ordered product of the effective theory operators describing the nonrelativistic 
QCD dynamics for the production and annihilation of the ti pair at leading logarithmic (LL) 
and NLL order. ^ At LL order (in dimensional regularization) the Green function has the 
simple analytic form 
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For the NLL order Green function we use the numerical techniques and codes of the TOPPIC 
program developed in Ref. 2J] (see also Ref. [25|) and determine an exact solution of the 
full NLL Schrodinger equation employing the approach of Refs. j21|. We estimate the QCD 
uncertainties in the normalization of the Higgs energy spectrum from the NLL order Green 



function as 5% [17|,|18|. Note that we account for the top quark finite lifetime by shifting 



the tt invariant mass Q used in the Green function into the complex plane such that the top 
quark relative velocity reads 



mt 



(8) 



where 



q2 = s + m]j-2^sE_ 



H 



(9) 



^ The renormalization scaling parameter v has mass dimension zero and is used in the effective theory to 
describe the correlated running of soft and ultrasoft fluctuations [l^l ■ The hard effective theory matching 
scale (at the top quark mass) is at i^ = 1 and low-energy matrix elements are evaluated for v r^ v ^ as to 
avoid the appearance of large logarithmic terms. 



This accounts for the top quark finite hfetime consistently at LL order, see for example |2f 
A consistent NLL description of finite lifetime effects and electroweak corrections shall be 
included in a subsequent publication. The term 5mt in Eq. (jH)) is a residual mass term that 
has to be specified perturbatively at each order to fix which top quark mass definition is 
being employed. In the pole mass scheme the residual mass term vanishes to all orders. We 
use the IS mass scheme J27|,|28|. The corresponding expression for 5mt at NLL order can also 



be found in Ref. jl7|. We use the IS top quark mass and implement the residual mass term 
in the soft factor of the factorization formula because it avoids the pole mass renormalon 
problem J29| and leads to a tt resonance peak position that is stable under higher order 
perturbative corrections J30|. For the NLL order QCD corrections to the hard factors, 
which are discussed below, we neglect the corrections that arise from the residual mass 
term because the numerical effects are at the 1% level and substantially smaller than the 
uncertainties from low-energy QCD effects. This approximation was also used in Ref. [17|. 
Concerning the hard contributions in Eq. (j2)), the first term in the parenthesis gives the 
contribution for the ti pair in a S-wave spin singlet state and the other three terms give the 
contributions for the tt pair in the three S-wave spin triplet (-1-1, 0, —1) states. As described 
already in Ref. jl7| we use the helicity basis for the top and antitop spinors in the endpoint 
where fci = k2 (see Fig. H]) to define the singlet and the triplet states. In this basis there 
are additional f-suppressed (NLL) contributions to the triplet ±1 contribution that arise 
from S-P wave interference terms and originate from the interference of vector and axial- 
vector contributions at the ti vertex. These additional order v contributions cancel in the 
sum of the triplet contributions and can also be avoided if a spin basis is used that does 
not depend on the momenta of the top quarks J3l|. Since here we are not interested in 
the phenomenology of top polarization these additional NLL order contributions are not 
included in Eq. (jS)). The functions F'^''^^ are the tree level (hard) matching conditions for 
the contributions of the respective ti spin states. They read 
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where 



Tl-2Qfsl Ti ^ e mt e e" ,^^, 

Here, Q/ and T^ are the fermion charge and weak isospin, e is the electric charge and s^ 
{cw) the sine (cosine) of the Weinberg angle. 

The functions Cj(z/) are the hard QCD Wilson coefficients and depend on mj, itlh and the 
cm. energy \/s. They also depend on the renormalization parameter v which accounts for 
the renormalization group running of the effective currents that produce and annihilate the 
ti pair in the various S-wave spin states. To achieve reliable predictions the renormalization 
scaling parameter u has to be chosen of order ctg, i.e. of order of the average top velocity in the 
tt cm. system. For this choice the imaginary part of the zero-distance Green function does 
not contain any large logarithms from ratios of the hard scales and the small nonrelativistic 
scales, the top three momentum p^ ~ mtv and the top kinetic energy Et ~ rritv'^ defined in 
the tt cm. system. All large logarithms are summed into the hard QCD coefficients. At NLL 
order the renormalization group evolution of the hard QCD coefficients can be parameterized 
as 

C[i,t),±{y) = C(i,,),±(l) exp(/(z/,2)) , (2 = 0,±1) 

co(z/) = co(l)exp(/(z/,0)). (14) 



The function / was given in Ref . [17[ using the results obtained in Refs. J20|, |32| • Whereas 
the renormalization group running of the coefficients can be determined within the effective 
theory, and is independent of the short distance process, the matching conditions at i/ = 1 
are process-dependent. We use the convention that the LL matching conditions for the 
Cj(z/) are normalized to unity. At NLL order the matching conditions are obtained from 
matching the factorization formula expanded to order as to the corresponding full theory 
Higgs energy distribution at C(as) in the endpoint region expanded to 0{y) for stable top 
quarks and using z/ = 1 (/i = mt) for the renormalization scaling parameters. The full theory 
predictions are taken from the numerical codes obtained in Ref. |15||. More information on 
the numerical matching procedure can be found in Ref. [l7|. The NLL matching conditions 
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TABLE I: Numerical values for the matching conditions for the singlet and triplet hard QCD 
coefficients for typical values ^/s, mt and mn- The masses and energies are given in units of GeV. 
Note that C(i +i) _|_ = C(i „i)^± due to parity. 



can be parameterized in the form 



ci_±(z/ = 1) = IH ()ci_±(Vs,mi,mH), 

Co(z/ = 1) = \^ dco{y/s,mt,mH), 



(^ = 0,±1) 



(15) 



and numerical results for the NLL order contributions for various choices of ^/s, mt and 
rriH are given in Tab. lU The singlet matching conditions do not depend on the electron- 
positron polarization because there is only one non-trivial QCD form factor in the full 
theory that can contribute to the hard QCD matching conditions for the effective theory 
spin singlet ti current. In Feynman gauge it originates from the pseudoscalar Goldstone-tt 
vertex. For the triplet currents, on the other hand, several form factors contribute in the 
full theory ti vertices, therefore the matching conditions are polarization-dependent for the 
parameterization used in Eq. Q. 
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FIG. 2: The unpolarized Higgs energy spectrum in the nonrelativistic expansion at LL (dotted 
hnes) and NLL (soUd hnes) order for v = 0.1, 0.2, 0.4. The fixed-order expansion is also shown at 
Born level (lower dotted line) and at 0{as) for /i = ^/s (lower dashed line) at for fi = ^/sv (upper 
dashed line). The cross section at NLL order fails to reproduce the correct physical behavior of 
the fixed-order results from the loop expansion in the low Higgs energy regime. At the IS peak 
the upper (lower) NLL order curve corresponds to the effective theory renormalization parameter 
i/ = 0.2(0.1). 



If the polarization of the tt final states is not accounted for, the factorization formula can 
be written in a simpler form using for the ti spin triplet contributions the definitions 



cu^)ni^ E 4.),±M^dl±' ni^ E ^d 
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cU-) = ^•^-'^'"■'b"''^'"^^''" . (16) 

In Ref. jl7| the results for the triplet contributions were presented in this form. 



III. THE LOW HIGGS ENERGY ENDPOINT REGION 

In Fig. 121 the prediction for the unpolarized Higgs energy spectrum obtained from the 
factorization formula in Eq. Q has been displayed at LL (dotted lines) and NLL (solid lines) 
order in the nonrelativistic expansion for the effective theory renormalization parameters 
1/ = 0.1, 0.2, 0.4. The parameters are ^i = 500 GeV, mf = 175 GeV, ttih = 120 GeV, and 



(17) 



We have also plotted the tree level (lower dotted line) and the 0{as) Higgs energy spectrum 
for /i = ^/s (lower dashed line) and for fi = ^v (upper dashed line) where v is the tt relative 



Tt = 1.43 GeV, 




Mz = 91.1876 GeV, 


Mw = 80.423 GeV 


«-! = 137.036, 


c^ = Mw/Mz . 



velocity defined in Eq. (jHJ. Since the hard scale as well as the relative momentum of the top 
quarks are scales that are relevant for nonrelativistic ti production, the difference between 
the two scale choices illustrates the ambiguity contained in the fixed-order calculation close 
to the large Higgs energy endpoint. A detailed discussion of the deficiencies of the fixed- 
order predictions in the endpoint region and quality of the nonrelativistic expansion and 
the theoretical normalization uncertainty of the NLL order prediction has been given in 
Ref. 113] and shall not be repeated here. The issue we want to point out in Fig. |21 is that 
the predictions obtained from the factorization formula in Eq. ©, which properly accounts 
for the summation of all NLL order contributions in the large Higgs energy region, is not 
compatible with the correct physical behavior at the low Higgs energy endpoint Eh = rnn- 
There the Higgs boson is produced at rest (in the lab frame) and the Higgs energy spectrum 
has to go to zero as do the tree level and 0{as) predictions. In particular, at the low Higgs 
energy endpoint region there is no singular enhancement from the matrix elements, and due 
to phase space suppression the coefficient functions Gi of e.g. the tree level Higgs energy 
spectrum (see Appendix^ vanish like Gi ^ /3 with 
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This endpoint behavior cannot be obtained within the nonrelativistic expansion in small v 
even if the endpoint is located at a velocity much smaller than one, see Eq. (j2)). Terms that 
are formally from beyond NLL order in v thus need to be summed up to achieve a correct 
low Higgs energy endpoint behavior. 

It is useful for the construction of a factorization formula which can account for the correct 
physical low Higgs energy behavior that the full theory tree level Higgs energy spectrum, 
both for the ti pair in the spin singlet and for the (combined) triplet configuration, does not 
have order v (NLL) corrections to the leading endpoint behavior in the large Higgs energy 
endpoint, i.e. 
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At NLL order it is thus consistent to use the full tree level Eh spectrum in the large Higgs 
energy endpoint instead of the constant LL matching conditions F^''^^ given in Eqs. ()1U|) - 






da 



n% 



dEHj^,,^Fi^^ + F2X 



2N,[{l+XH-Axtf-AxH, 

_g3/2^ 



1/2 



0,± 



da 



T?Z 



dEHJs,,^F,^^^ + Fll 



2N, [{l+XH-4:Xtf-4:XH 

gZ/2 ^ 



1/2 



(20) 



where {-^-)^orn ^^ ^^^ ^^^^ ^^^^ level Higgs energy spectrum for the polarized e~^e initial 



state. Note that the replacement prescription in Eq. (|^ can only be applied for Higgs 
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energies smaller than E^, for larger Higgs energies Eq. Q is left unchanged. For the 
convenience of the reader we have given the analytic expressions for the tree level Higgs 
energy spectrum in the appendix using up to minor modifications the conventions of Ref . |13| . 
They also correct a few typos that were contained in Ref. [ISf and pointed out before in 
Ref. [17[. For the case of an unpolarized tt final state, using the prescription (^01) in the 
factorization formula (jH)) leads to a modified factorization formula that resums correctly all 
NLL order terms. In addition it has the correct physical behavior at the low Higgs energy 
endpoint Eh = rnu- 

The modified NLL factorization formula based on Eqs. Q, (fTBj) . and (PUJ) is not unique, 
alternative prescriptions to achieve the correct physical low Higgs energy endpoint behav- 
ior are conceivable. However, different prescriptions will only affect the low Higgs energy 
endpoint where the Eh spectrum vanishes, and they should therefore not have a large nu- 
merical impact. While the modified NLL factorization formula contains the exact tree level 
contribution, its 0{as) contribution (in the expansion in powers of a^) differs from the exact 
0{as) result obtained in Ref. [15| since it includes only the QCD corrections of the large 
Higgs energy endpoint. Thus an estimate of the intrinsic uncertainty in our prescription 
can be gained by comparing its 0{as) terms with the exact result from Refs. [Il|,|l5|. For 
stable and unpolarized top quarks the first two terms in the as expansion of our modified 
factorization formula read 



da 



dE 



H 



{Eh, 



NLL 



da 



dE 



H 



{Eh, 



+ 



Born 



da 



dE 



H 



{Eh, 



+ 0{a 



S> 5 



OK) 



(21) 



where 



da 



dEn 



{Eh, 



± 



0[as) 



Cpas 



EL + E^' 



4m? 
1 



-1/2 



da 



dE 



{Eh, 



H 



Born 



(22) 



In Tab. ITTlnumerical results for the exact total 0{as) unpolarized cross section, a, 



exact 



and for the 0{as) approximation from Eq. (PT|) . a^i^\ are shown for various cm. energies 
and rut = 175 GeV, niH = 120 GeV, Fj = 0, /i = ^/s. For cm. energies below 500 GeV 
the deviation increases with the cm. energy. It vanishes at the three-body threshold ^/s ~ 
2mt + rriH and reaches the level of 1.5% for ^/s = 500 GeV. 

In Fig.inithe exact 0{as) unpolarized Higgs energy spectrum (black lines) and the 0{as) 
approximation in Eq. (PT| (gray lines) are displayed in 0.1 GeV bins for ^/s = 490, 500, 600, 
and 700 GeV, rrit = 175 GeV, iriH = 120 GeV, Ft = 0, and /i = a/s. Note that for the strong 
coupling we use as(500 GeV) = 0.09396. The other parameters are chosen as in Eq. (jTTj). 
For ^/s = 500 GeV the relative deviation in the Higgs energy spectrum is at most 2.8%. The 
difference is smaller for lower cm. energies since the maximal possible top relative velocity 
^max jg increasing with the cm. energy, see Eq. Q. The results indicate that the intrinsic 
uncertainty of our approach is substantially smaller than the theoretical uncertainty of 5% 
from uncalculated higher order QCD effects |l7l.ll8|. 

In Figs.Elc and d and in Tab.lTTlwe have analyzed the difference between the exact 0{as) 
results and the 0{as) approximation based on Eq. ()2H1 for larger cm. energies as well. It is 
a surprising fact that the fairly simple expression in Eq. (J21I), which contains only tree level 
information and the NLL QCD information from the large Higgs energy endpoint, can also 
account very well for the exact 0{as) results at higher energies, where real gluon radiation 
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V~s 


'^ exact 


^NLL 


rel. dev. (%) 


475 


0.0311 


0.0309 


0.6 


480 


0.0908 


0.0900 


0.9 


490 


0.254 


0.251 


1.2 


500 


0.446 


0.439 


1.5 


550 


1.366 


1.343 


1.7 


600 


1.953 


1.924 


1.5 


700 


2.356 


2.348 


0.4 



TABLE II: The total cross section using the exact 0{as) result cr^^^^l from Ref. [15| and the 
approximation based on Eq. H21() . Cj^ll" • '^^^ third column shows the relative deviation in percent. 
The difference between the two calculations is maximal for cm. energies around 550 GeV. 



/s = 490 GcV 




119 121 123 125 127 

E„ (GeV) 
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v^ = 
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0.030 
0.025 
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- 


—■■ 0.015 
^ 0.010 
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-_ 
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FIG. 3: The exact 0{as) unpolarized Higgs energy spectrum from Ref. [l5| (black lines) and the 
0{as) approximation in Eq. (PT|l (gray lines) for different cm. energies y/s for nit = 175 GeV, 
niH = 120 GeV, and /x = ^/s. 



is non-negligible. For cm. energies between 500 and 700 GeV the approximation based on 
Eq. (|7Hl deviates from the exact results by at most 1.8% for the unpolarized total cross 
section, where the maximal deviation is reached for ^/s ^ 550 GeV. Since the numerical 
evaluation of Eq. (j^ is substantially faster than for the exact 0{as) result [15[, it can be 
useful as an efficient approximation formula for higher cm. energies. 
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; m,"^ = l7SGeV. m„ = ]20GeV 
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Eh (GeV) 




m,'^=175GeV, m„=]20GeV 
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3, 0.01 
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m!^ = l75GeV, m„=J20GeV_ 
^/l=500GeV^^. 




Eh (GeV) 
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FIG. 4: The unpolarized Higgs energy spectrum for different cm. energies at NLL order (solid lines) 
using the modified factorization formula based on Eqs. Q, (|16|) . and (|20|) for the renormalization 
parameters u = 0.1,0.2,0.4, at 0{as) (dashed lines) from Ref. [uj with fi = ^/s, and at Born 
level (dotted line). At the IS peak the upper (lower) NLL order curve corresponds to the effective 
theory renormalization parameter u = 0.2(0.1). 



IV. NUMERICAL ANALYSIS 



In Fig. m the unpolarized Higgs energy spectrum at NLL order (solid lines) using the 
modified factorization formula based on Eqs. Q, ()16|). and ()20j) is displayed for the renor- 
malization parameters u = 0.1, 0.2, 0.4 for the cm. energies y/s = 485, 490, 495, 500 GeV and 
rrit = mj^ = 175 GeV, rriH = 120 GeV. The other parameters are chosen as in Eq. (|17|). For 
comparison we also show the tree level prediction (dotted lines) and the 0{as) results [15| 
(dashed lines) with /i = y^ for a stable top quark. The nonrelativistic NLL order results 
show a substantial enhancement compared to the tree level and one-loop QCD predictions. 
The Higgs energy spectrum in the effective theory extends beyond the endpoint E^ that 
is obtained for the stable top quark case. This is because the top quarks can be produced 
off-shell with invariant masses smaller than rrit if the top quark decay is accounted for. 
With the present technology the finite top quark lifetime can only be implemented system- 
atically in an expansion in the top quark off-shellness, which is naturally provided by the 
nonrelativistic expansion we use here. 

It is conspicuous that the spectrum above the endpoint E^ in the NLL prediction falls off 
quite slowly. Since the average cm. top quark velocity increases with the Higgs energy for 
Eh > E'jj we define the total cross section by applying a cut on the Higgs energy above E^ 
such that the average cm. top velocity remains below t;cut = 0.2. We fix the relation between 
the maximal Higgs energy and Wcut by the relation E'^* = (s + m|^ — Ql^^)/{2y/s), which is 
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v^ 


mt 


niH 


^^orn(fb) 


4LL(fb) 


'^NLL/'^Born 


^Born(fb) 


^NLL(fb) 


'^NLL/'^Born 


500 


170 


115 


0.644 


0.989(49) 


1.54 


1.660 


2.568(128) 


1.55 


490 


170 


115 


0.444 


0.754(37) 


1.70 


1.149 


1.965(98) 


1.71 


480 


170 


115 


0.260 


0.516(25) 


1.98 


0.674 


1.347(67) 


2.00 


470 


170 


115 


0.108 


0.285(14) 


2.64 


0.281 


0.747(37) 


2.66 


460 


170 


115 


0.014 


0.086(4) 


6.17 


0.036 


0.226(11) 


6.21 


500 


170 


120 


0.486 


0.783(39) 


1.61 


1.258 


2.040(101) 


1.62 


490 


170 


120 


0.312 


0.568(28) 


1.82 


0.809 


1.483(74) 


1.83 


480 


170 


120 


0.159 


0.355(17) 


2.23 


0.413 


0.929(46) 


2.25 


470 


170 


120 


0.046 


0.159(7) 


3.48 


0.120 


0.418(20) 


3.50 


500 


170 


140 


0.102 


0.229(11) 


2.24 


0.268 


0.604(30) 


2.26 


490 


170 


140 


0.029 


0.101(5) 


3.48 


0.076 


0.268(13) 


3.51 


500 


175 


115 


0.459 


0.787(39) 


1.72 


1.181 


2.039(101) 


1.73 


490 


175 


115 


0.268 


0.538(26) 


2.01 


0.692 


1.399(69) 


2.02 


480 


175 


115 


0.111 


0.298(14) 


2.68 


0.288 


0.777(38) 


2.70 


470 


175 


115 


0.014 


0.091(4) 


6.32 


0.037 


0.236(11) 


6.35 


500 


175 


120 


0.322 


0.593(29) 


1.84 


0.832 


1.541(77) 


1.85 


490 


175 


120 


0.164 


0.371(18) 


2.26 


0.425 


0.967(48) 


2.28 


480 


175 


120 


0.047 


0.167(8) 


3.54 


0.123 


0.437(21) 


3.56 


500 


175 


140 


0.030 


0.107(5) 


3.55 


0.079 


0.281(14) 


3.57 



TABLE III: The total cross section in units of fb at Born level for stable top quarks and at NLL 
order for unstable top quarks using z^ = 0.2 for fully polarized electron-positron beams. The index 
refers to the polarization of the electron beam. The masses and ^/s are given in units of GeV. For 
mt = (170, 175) GeV we use Tt = (1.31, 1.43) GeV. 



exact in the stable top case. Here, Ql^^ = (4mj)/(l + fcut) is the minimal tt invariant mass. 
Note that Qcut is smaller than 2mt because for E^ > E^ we are in the bound state regime. 
As mentioned before, we plan a systematic treatment of finite lifetime and off-shell effects 
at the NLL order level in a subsequent publication. 

In Tab. Illll the impact of the NLL order summations on the total cross section for un- 
polarized ti pairs and polarized electron-positron beams is analyzed numerically for various 
cm. energies, top quark masses and Higgs masses. The other parameters are chosen as in 
Eq. (fT7|) except for the case mt = 170 GeV where we use Tt = 1.31 GeV. In Tab. IllH 
<7Born refers to the tree level cross section for stable top quarks (see the appendix for explicit 
expressions) and ctnll to the NLL total cross section as defined above and based on the 
modified factorization formula discussed in Sec. IIIII The NLL order predictions were ob- 
tained for the effective theory renormalization parameter u = 0.2. The uncertainties given 
for (Tnll reflect the 5% theoretical error from higher order QCD and relativistic corrections 
as discussed in Ref. jl7|. The results in Tab. IIIII demonstrate the importance of the summa- 
tion of the singular terms proportional to (Os/f)" and (a^lnf)" that arise in the endpoint 
region, and of the off-shell effects that arise from the finite top quark lifetime. Compared to 
the tree level predictions the enhancement is more pronounced for smaller cm. energies and 
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FIG. 5: The total cross section for unpolarized top quarks at tree level (dashed lines) and at 
NLL order (solid lines) as a function of ^/s (left panel) and as a function of rriH (right panel) for 
unpolarized electron-positron beams (P+, P-) = (0, 0) (respective lower curves) and for (P+, P-) = 
(0.6, —0.8) (respective upper curves). 



larger top or Higgs masses. 

It is a realistic option for the ILC project to polarize the e~^e~ beams up to {Pj^.,P-) = 
(0.6,-0.8) p. Since this can further enhance the cross section we have also assessed its 
merits for the process at hand. In Figs. the total cross section for unpolarized top quarks 
at the tree level (dashed lines) and at NLL order (solid lines) is shown as a function of y/s and 
rriH for unpolarized electron-positron beams (P+, P_) = (0, 0) and for (P+, P_) = (0.6, —0.8). 
The other parameters are chosen as in Eq. (fTTj) . see also the figure caption for more details. 
For the NLL cross section the predictions for the three choices z/ = 0.1,0.2,0.4 for the 
renormalization scaling parameter are shown. The results demonstrate that using electron- 
positron polarization the cross section can be enhanced by roughly a factor of two over the 
unpolarized cross section. Compared to the tree level predictions for unpolarized electron- 
positron beams, which were the basis of previous experimental analyses '22], QCD effects 
and beam polarization {P^,P_) = (0.6,-0.8) can enhance the cross section by about a 
factor of 4 or even more for ^/s = 500 GeV, depending on the Higgs mass. Because of the 
limited statistics for tiH production during the first phase of the ILC project, these results 
are important for realistic experimental simulations of Yukawa coupling measurements. 



V. CONCLUSION 



We have analyzed the impact of summing the QCD singularities proportional to (cts/f )" 
and {as Inf )" that arise in the large Higgs energy endpoint region for the process e+e~ — > tiH 
for cm. energies up to 500 GeV, i.e. energies which can be achieved during the first phase 
of the ILC project. The singularities cause the breakdown of usual multi-loop perturba- 
tion theory in powers of Og and originate from nonrelativistic dynamical QCD effects that 
arise because the relative velocity of the ti pair is small. A consistent theoretical treat- 
ment requires the use of nonrelativistic effective theory methods and includes a systematic 
treatment of off-shell effects caused by the finite top quark lifetime. In Ref. |17| we derived 
a factorization formula for the large Higgs energy endpoint region for large cm. energies 
above 500 GeV. In the present work we have extended the approach to cm. energies below 
500 GeV, where the top quark pair is nonrelativistic in the entire phase space, and we have 
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also accounted for the effects of electron-positron beam polarization. We have determined 
the predictions for the Higgs energy spectrum and the total cross section at NLL order for 
the QCD effects and at LL order for the top quark finite lifetime and for off-shell effects. The 
QCD effects enhance the total cross section by roughly a factor of two relative to the Born 
prediction for ^/s = 500 GeV. Using polarized electron-positron beams the cross section can 
be further enhanced over the unpolarized case by another factor of approximately two. Our 
results are important for realistic simulation studies for Yukawa coupling measurements in 
the first phase of the ILC project. 
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APPENDIX A: TREE LEVEL HIGGS ENERGY SPECTRUM 



Correcting the typos of Ref. jl3| the tree level Higgs energy spectrum in the process 

e~^e~ — > tiH for polarized electron-positron beams reads {xe = '^^En/y/s, Upt = 47ra;^/(3s)) 



fdaiEnlY _ ^ 
V dxE J.„_„ " """"'Sn^ 



Born 



1-Xz 



QlQ\ + 



+ 



+ 



[Ve =F QeJ 



^'J2G^ + vKG, + Ge 



i=2 



+ 



{l-xzY 

QeQtiVeTae)Vt 
1-Xz 



Gi 



G 



6 f 5 

(Al) 



where the coefficient functions are given by 



Gi 



2A? 



,^2 2N ■ -^K^^t-XH){2Xt + l)XE 

[P ~ X^jXe 

+ (/^^ - x\) [iQxl + 2x]j - 2xhXe + x\- 4xt(3xH - 2 - 2xe)\ In ( -^ 

\xe - 



Xe + P 
xe - (3 



G, 



-2A? 



f3xE [ — 96x^ + 24:XtXH — {—xh + 1 + xe){x'^e ~ P' 



(/32 - x%)xe 

+ 20^ - x|) [24x^ + 2{xl - xhXe) + Xt{-lAxH + 12xe + x|)] In ( ^^^ 

\xe - 



Xe+ P_ 



(A2) 



(A3) 



These first two coefficients describe the s-channel exchange of the photon and the Z boson 
where the Higgs boson is radiated off one of the top quarks |l3|. A missing factor s is 
introduced in the first line of the formula for G2- 
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The coefficient functions G3 to Gq describe the emission of the Higgs boson from the 
Z-boson, 

G3 = —, ~'^^^f' 1 4x1 + 12x| + 2xzxl 

xz{xh - xz + 1- xeY y 

+ {-1+xe)xI-xh[^xz + {-4. + Axe + xI)\ I, (A4) 

^^ = ^i ^^^^1 y^ (48^* + 12a;^ - (-24 + /^^ + 240;^, - 3x|) j , (A5) 

G5 = — ;^ ^-^-^ II3[Qxz + xe{-xh-^ + xe)\ 

Xz {-^H +Xz -1+xe) I 

+ 2 [x^(x^ - 3x2 + 1 - xe) + Xt{-4xH + 12x2 + 4)] In ( ^^^ ) | , (A6) 

\xe -13/ J 

G. = -^>-fi->-f'' h + (4., - .„ + 2 - X,) In (££±4) I . (A7) 

-xh+xz-I + xeI \xe-(3J) 

These terms give contributions to the Higgs energy spectrum of less than a few percent in 
the energy range between 500 GeV and 1 TeV. The overall signs of G5 and Gg are changed 
relative to [13J. The couplings and constants are defined in Eqs. ()13I4|) and the term /3 is 
given by 

^ \V~siiE%-EH)yrs + 2m^)J ' ^ *^ 

with the large Higgs energy endpoint being defined as 

_ s + mj,- Ami .,Q^ 

En = 1^^ ■ (A9) 
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